arXiv:1502.07480v2 [cond-mat.mtrl-sci] 11 Ju12015 
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We find that the motion of the valley electrons - electronic states close to the K and K' points of 
the Brillouin zone - is confined into two dimensions when the layers of M0S2 form the 3 R stacking, 
while in the 2 H polytype the bands have dispersion in all the three dimensions. According to 
our first-principles band structure calculations, the valley states have no interlayer hopping in 3 R- 
M0S2, which is proved to be the consequence of the rotational symmetry of the Bloch functions. By 
measuring the reflectivity spectra and analyzing an anisotropic hydrogen atom model, we confirm 
that the valley excitons in 3R-M0S2 have two-dimensional hydrogen-like spectral series, and the 
spreads of the wave function are smaller than the interlayer distance. In contrast, the valley excitons 
in 2H-M0S2 are well described by the three-dimensional model and thus not confined in a single 
layer. Our results indicate that the dimensionality of the valley degree of freedom can be controlled 
simply by the stacking geometry, which can be utilized in future valleytronics. 

PACS numbers: 73.22.-f, 71.20.Nr, 71.35.-y, 78.20.-e 


I. INTRODUCTION 

The monolayer systems of group-VI dichalcogenides, 
MX 2 (M=Mo, W; X=S, Se, Te) [I] have received con¬ 
siderable interest as unique alternatives of graphene 
for their various intriguing properties [2]: direct band 
gap of 2 eV [3HT2], high photoluminescence yield [31 
0 ], high on/off switching ratio in field effect transis¬ 
tors mm , and electric field-induced superconductiv¬ 
ity 05] 03]. Most notably, their staggered-honeycomb- 
like lattice structure (the structure and the unit cell 
are shown in Fig. [l] (a)-(c)) hosts a model system for 
valley-dependent phenomena m originally proposed for 
graphene [731 - 1301 . The bottom of the conduction band 
and the top of the valence band are located at the K 
points of the hexagonal Brillouin zone in MX 2 . Since 
the K and K'= —K points are not equivalent, the elec¬ 
tronic states in the opposite pockets can carry an addi¬ 
tional quantum number, the valley index. Furthermore, 
the spin-orbit coupling results in sizable valley-dependent 
spin splitting at the valence top (VT). On the basis of 
the strong coupling between the spin, orbital and valley 
degrees of freedom, the control of the carrier population 
of each valley by orbital/spin-sensitive probes has been 
proposed. Indeed, valley-selective excitation of the elec¬ 
trons and excitons by circularly polarized light has been 
demonstrated [27137)1 . Further understanding and sta¬ 
ble control of the valley-dependent electronic properties 
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could bring us to the valleytronics —a way of encoding 
information into the valley degree of freedom, which can 
be much faster and more efficient than conventional op¬ 
toelectronics H7. [781126] . 

The above progress led to renewed attention to the rich 
polymorphism of multilayered MX 2 [3 [27] in view of the 
valley physics. In inversion-symmetric bilayers with the 
2H stacking [Fig. [I] (d)], the net valley-dependent spin 
polarization is absent. This property has been proposed 
to be utilizable for switching of the polarization with a 
symmetry-breaking electric field |28l 23] or mechanical 
bending [S3]. It has also been stated that a layer degree 
of freedom (upper/lower layer) couples to the valley de¬ 
grees of freedom mm and this coupling should cause 
magnetoelectric effects m- On the other hand, very re¬ 
cently, Suzuki et al. [33] observed valley-dependent spin 
polarization in multilayered M 0 S 2 by utilizing the non¬ 
centrosymmetric 3R stacking [Fig. [T| (e)]. This success 
paves a very different way to the control of the valley car¬ 
riers: Valleytronics by engineering the stacking geome¬ 
try. However, knowledge of the valley electronic states in 
the 3R systems is still scarce, which hampers further ad¬ 
vances in the field and applications of multilayered M 0 S 2 . 

In this article, we study the valley electronic states 
in M 0 S 2 with the 3R stacking and compare them with 
those in 2H-MoS2- Combining ab initio band structure 
calculations and group-theoretical analysis, we show that 
the interlayer hopping amplitude of the valley states is 
exactly zero in 3 R-M 0 S 2 , i.e. the electrons are confined 
within the two-dimensional (2D) layers. Furthermore, we 
study how this confinement affects the exciton spectrum 
with an anisotropic hydrogen atom model. Finally, we 
compare the theoretical results to the reflectivity spectra 
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FIG. 1 : (a) Trigonal prismatic unit cell of monolayer M0S2; 
(b)-(c) side and top views of the monolayer, where the trig¬ 
onal prisms are depicted as shaded triangles, a, /3 and 7 
represent inequivalent three-fold rotational axes and ai and 
a2 are primitive lattice vectors, (d)-(e) side view (top) and 
top view (bottom) of the 2 H and 3 R stackings, respectively. 
Rotational axes for each layer are also indicated. 




FIG. 2: (a)-(b) Band structures of bulk M0S2 with the 

2 H and 3 R stacking, respectively, where we employed spe¬ 
cial points of the conventional hexagonal Brillouin zone (BZ) 
defined by r=(0,0,0), A=( 0 , 0 ,7r/c), K=( 0 ,47r/(3a), 0 ) and 
H=( 0 ,47r/(3a),7r/c) with a=|ai|=|a2| and c being the inter¬ 
layer distance. Note that the the present BZ is twice as large 
as the primitive BZ of 2H-M0S2. Solid (dashed) lines denote 
the result with (without) spin-orbit coupling. 


as implemented in the WIEN2k code. Experimental lat¬ 
tice parameters and atomic configurations were taken 
from Refs. [551 and (55). The muffin-tin radii for Mo and S 
atoms, tmo and rs, were set to 2.44 and 2.10a.u., respec¬ 
tively. The maximum modulus for the reciprocal lattice 
vectors AT max was chosen so that rsxA' max =7.00. 

The calculated band structures are shown in Fig. [2] 
The apparent difference seen around the T point in va¬ 
lence bands is due to Brillouin-zone folding for the 2H 
case with respect to fc z =7r/(2c) plane and hence trivial. 
Notably, along the K H path, the band at the conduction 
bottom (CB) is flat for the both polytypes. On the other 
hand, the VT bands for the 2H polytype shows small but 
nonzero dispersion, whereas those for 3R is flat. This 
feature has been first found in Ref. 133] but its origin has 
remained unexplored. Below, we clarify its mechanism 
analytically. 


measured on both 3 R-M 0 S 2 and 2 H-M 0 S 2 compounds. 
The revealed mechanism of the single-layer confinement 
of the valley electrons respects only the crystal and or¬ 
bital symmetries and therefore is commonly applicable 
to the family of 3 R-MX 2 , which should facilitate the di¬ 
mensionality control of the valley states in the transition- 
metal dichalcogenides. 


II. THEORY ON THE VALLEY ELECTRONIC 
STATES 

A. First-principles calculation 

First, we calculated the band structures of bulk 
2H- and 3 R-M 0 S 2 using the WIEN2k code employ¬ 
ing the full-potential linearized augmented plane-wave 
method [34] . We used the Perdew-Burke-Ernzerhof 
exchange-correlation functional [35] and included the 
scalar-relativistic effects [36] and spin-orbit coupling EH 


B. Group-theoretical analysis 

In order to understand the microscopic origin of the 
flat bands along the K-H path, let us first look into the 
symmetries of the electronic states in a monolayer of type 
A in Fig. |T] The little group of K point includes three fold 
rotation (C 3 ). Under the rotation around the axis passing 
through the center of the honeycomb unit [a in Fig. Boi, 
the VT (CB) state, which is formed by Mo Ad x 2 _ y 2 +d xy 
( 4 d 2 2 ) orbital having orbital angular momentum l z =—2 
( 0 ), gets phase 1 (e*^) [?51 . 

As a matter of fact, there are two more rotational axes 
inequivalent to a: /? and 7 passing through S and Mo 
sites, respectively [Fig. |TJ)c)] . In order to calculate the 
phase factors for these C 3 axes, let us decompose the C 3 
rotation into two successive operations: permutation of 
the atomic sites P and rotation of the local coordinates 
at each lattice point R. The C 3 around j axis is defined 
as C 3 .j = RPj (j = a,/?, 7), where P } denotes the corre¬ 
sponding site permutation. The contribution to the total 
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FIG. 3: (a)-(c) Schematic of the phase change due to site 
permutation. Real-space configuration of the phase of the K- 
point Bloch state [exp(iK ■ r)] is represented, where circle, 

- 27T 

triangle and square represent the three phase values 1 , e J 3 

• 47T 

and e l 3 . respectively. Large shaded triangles correspond to 
those in Fig [T] Phase configurations before (after) the cor¬ 
responding permutation Pj are represented by solid (shaded) 
symbols. Crosses indicate the rotational centers j = a, p and 
7 , whereas arrows represent how the sites permute with Pj. 


phase from R is simply with 0=2n/3. The contri¬ 
bution from Pj can be calculated by examining how the 
Bloch wave function transforms under the site permu¬ 
tation. In Fig. [3j the phase of the K-jooint Bloch state 
is depicted, where three values ( 1 , e 1 ^ and e'f) are 
represented by circle, triangle and square, respectively. 
By operating P a , e.g., the configuration depicted with 
solid symbols is transformed to that with shaded sym¬ 
bols [panel (a)]. Comparing these symbols on each site, 
the phase contribution e l ~z for P a is deduced. Applying 
similar arguments to j=/3 and 7 , and to the monolayer of 
type A\ we get the results summarized in Table [I] Note 
that the phases for A are inverse of those for A because 
their valley states are related by spatial inversion and 
time reversal. 

On the basis of the axis-dependent phases derived 
above [40], let us next consider the interlayer hopping 
amplitude defined by 

t b =(^ b KA \n\^ 2 ) (b = VT, CB) ( 1 ) 

for the bulk material. Here, |'F^ C L ) denotes the K-point 
Bloch functions at Lth layer, which are eigenstates in 
the limit where the layer is isolated from others. R is 
the total Hamiltonian of the bulk crystal. Suppose C 3 
denotes a certain three-fold rotation under which R is 


TABLE I: Phases Sj of the eigenvalues e 2nlS j of C 37 = RPj 
around different rotational axes (j =a, p and 7 ) for the K- 
point states ( 6 =VT, CB). 
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invariant, then 

(^K,ll^l^K,2) = (^, 1 l C 3“ 1 C3^C 3 - 1 C3|^ i2 ) 
=(«i|C' 3 - 1 )^(C 3;/ |^, 2 )) 

=ex p[—27r*Ai 2; jj'](^ , K,il’^l^ r K,2)- (2) 

Here, - ^ 27 ' an d denotes the phase sum¬ 

marized in Table[l]for the Lth layer. Note that j and j' in 
the second and last lines depend on the stacking pattern 
and rotational axis of C 3 . In the 2H stacking, without loss 
of generality, we can assume the first and second layers to 
be A and A', respectively. In this case, the possible com¬ 
binations (j,f) are {a, a), (f3, 7 ) and (7 ,/3) [Fig. 0 (d)], 
all of which yield A 22 .^.,=0 ( 6 =VT) and 2/3 (&=CB), 
respectively. In the 3R stacking, where only layers of A 
type are stacked, (j,j')=(a, p), (/3, 7 ), (7 ,a) [Fig. 0(e)], 
and therefore Aj 2 ; ^-,=l/3 (b = VT, CB). According to 
Eq. pi, the nonzero AJ 2 ■ forces R to be zero, which 
consistently explains the dispersionless bands along the 
K-H direction in the 3R polytype. 


III. VALLEY EXCITONS 
A. Experiment 


The difference in the band structures of 2H- and 3R- 
M 0 S 2 should cause different optical properties of their 
valley excitations. To verify this experimentally, we mea¬ 
sured the normal incidence reflectivity spectra on the 
cleaved (001) surfaces of both polytypes at T=20K (see 
Fig- 0)- The details of the crystal growth were reported 
in Ref. [-TT . Since the reflectivity were measured in the 
limited photon energy range of £’=1.5-3 eV, the spectra 
were extended up to 30 eV using the data published for 
2 H-M 0 S 2 [43] to facilitate proper Kramers-Kronig trans¬ 
formation. 

The imaginary part of the dielectric function calculated 
from the reflectivity is shown in Fig. 0 Although the 
two main features in the £2 spectra, the A and B exciton 
peaks corresponding to the optical excitations from the 
spin-orbit split VT to the CB, look similar, there are 
several important differences between the two spectra. 
First, the energy separation between the A\ and the B 
peaks is larger for the 2H compound (197 meV) compared 
with the 3R polytype (148meV), which agrees with the 
theoretical prediction for the splitting at the VT (see 
Fig. 0and Ref. 133]) . Beside the main excitonic resonances, 
smaller peaks were observed on the high energy shoulder 
of the A exciton, labeled as A 2 . In a hydrogen-atom-like 
model of excitons, A\ and A 2 are assigned to the Is and 
2s bound states of the electron-hole pair. In the 3R case, 
Ai and A 2 peaks are more separated and their intensity 
ratio is larger than in the 2H structure. 
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FIG. 4: (a)-(b) Reflectivity spectra of bulk 2H- and 3R- 
M 0 S 2 measured at T=20K. The imaginary part of the di¬ 
electric function was obtained from the reflectivity spectra by 
Kramers-Kronig transformation. For the A exciton branch, 
Is and 2s excitons are labeled as Ai and A 2 , while the fine 
structure of the B exciton cannot be resolved. 


B. Model analysis 

To understand these differences, we employed an 
anisotropic hydrogen atom model 1111 : 


v£+vg 

2to* 


V|_ 

2tom 


1 




[x 2 + y 2 )+e 2 ± 

= Ecp( r), 


V?(r) 


(3) 


where £± and £y are the components of the dielectric con¬ 
stant tensor perpendicular and parallel to the ^-direction 
respectively, and m* j_ and my are those of the reduced 
effective-mass tensor of excitons at the K-point. Note 
that in a monolayer system, the screened electron-hole 
interaction is rather different from that in Eq. (|3| ow¬ 
ing to the surrounding vacuum |451146| . Bulk environ¬ 
ment allows us to use the present model as in early stud¬ 
ies EH SB]. 

TABLE II: Parameters adopted for Eq. ([ 3 ] extracted from ab 
initio band structure calculations: dielectric constants and 
reduced effective masses of excitons at the K-point (in units 
of electron mass mo). 


The values of e±, £y, m ^ and m| were extracted from 
ab initio calculations using the WIEN2k code. Dielectric 
constants £ were calculated using the random-phase ap¬ 
proximation 


c 1 occ.unocc.i / , 

, = i + | Hm iyy 


% 


'zkcr / 


^■ikcr 


(4) 


where V denotes the volume of the simulation cell, i and 
j are indices of the occupied and unoccupied bands, re¬ 
spectively, k is the wave vector, e a (a = x, y , z) denotes 
the unit vectors, a denotes the spin index, fa k CT denotes 
the Kohn-Sham wave functions, and £jko- denotes their 
orbital energies. The effective masses m* were estimated 
from a parabolic curve fitting of the band dispersions of 
electrons and holes around the K point within the range 
of ~ 0.02 A -1 (~ 0.03 A- 1 for my of the 3R polytype). 
To obtain the effective masses accurately, r m 0) t s, and 
rsxKT max were set to 2.20a.u., 1.89a.u., and 9.00, re¬ 
spectively. 

Calculated dielectric constants and reduced effective 
masses are presented in Table JTT] which agree well with 
previous calculations [52HM]- Since both VT and CB 
are flat along the K H line, Wy is quite large in the 3R 
stacking, whereas it is small in the 2H stacking in which 
only CB is flat. Note that the effective masses were de¬ 
termined in the close vicinity of the K point, and the va¬ 
lence and conduction bands in 3 R-M 0 S 2 have less than 
2 meV dispersion along the K-H line. Since this range is 
much smaller than the binding energy of Is or 2s exci¬ 
ton, in the following analyses, we regard the 3R bands 
completely flat, i.e. TOy=oo. 

Within the model of Eq. ([ 3 ]), the binding energies and 
oscillator strengths of the Is and 2s excitons were esti¬ 
mated. For the 2H case, we referred to a previous report 
presenting a numerical solution of the model for a range 
of parameters [48]. For the 3R case, we instead took 
the exact solution of the 2D hydrogen atom according 
to the above consideration. The ratio of the oscillator 
strengths was derived with the Fermi golden rule |55j 
using the present energy eigenvalues and the reduced 


TABLE III: Exciton binding energies of Is exciton Ei 3 , 2s ex¬ 
citon E 2s , their difference A E = E 2s — E\ a (all in meV), and 
a relative ratio of oscillator strengths of Is and 2s excitons 
fis/fis- Theoretical values are obtained from the anisotropic 
hydrogen atom model (Eq. <§)• Experimental values are ob¬ 
tained from a Lorentz model fit for the reflectivity spectra 
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effective masses. The obtained values are presented in 
Table |TTT| together with our experimental ones extracted 
from fits with Lorentzian peaks (Fig. [4]). The substan¬ 
tial differences in the energy levels and in the oscillator 
strengths between the exciton spectra of the 2H and 3R 
structures were captured by the model calculation. Since 
all the parameters in the 2H and 3R cases are similar to 
each other except m|, we attribute the differences solely 
to the 2D confinement of the electrons and holes in the 
3R compound. 

To see another aspect of the 2H and 3R excitons with 
respect to their dimensionality, we examined the exciton 
radii by using a simple trial wave function [56[ for the Is 
eigenstate of Eq. 


<Pis(r) 



-sj{x 2 + y 2 )/r]_ + z 2 / r fj 


(5) 


with the parameters r_L and ry, which represent exciton 
in-plane and out-of-plane radii, respectively, optimized in 
a variational manner. For the 3R exciton, we used a finite 
mj| listed in Table |Tlj to estimate an upper bound of ry. 
The obtained values of rjj_ are almost the same values for 
the 2H and 3R structures: 19 and 15 A, respectively. In 
contrast, the values of ry are 8.9 and 2.7 A for the 2H and 
3R structures, respectively. Comparing these values with 
the interlayer distance (6.1 A for both structures), we can 
see that the 3R exciton is tightly confined in a single layer 
while the 2H exciton extends to the neighboring layers. 


have much longer lifetime in 3 R-M 0 S 2 . Now we have 
clarified that the interlayer hopping amplitude is exactly 
zero and thus the protection is robust; also, it is com¬ 
mon to all 3 R-MN 2 . Note that the 2D states in the 3R 
systems are different from those in the monolayer sys¬ 
tems. In the former, the thickness of the sample allows 
us to generate valley polarized carriers away from the 
substrate, which will reduce impacts of the substrate po¬ 
tential. Hence, for example, we can use the 3R crystals 
as a medium (or lead wires, e.g.) for the valley device, 
through which the valley-polarized carriers are robustly 
transmitted. The unique property of the 3R stacking 
could thus bring about more effective valley-information 
transmission and generation toward dissipation-less low- 
energy optoelectronic devices. 


V. SUMMARY 

To summarize, we showed rigorous absence of the inter¬ 
layer hopping for the valley electrons in 3R-M0S2, which 
was substantiated with a group theoretical argument. 
Through the measurement of the reflectivity spectra and 
analysis with an anisotropic hydrogen atom model, the 
absence of the hopping has been shown to make the exci¬ 
ton absorption spectrum 2D hydrogen-like and the exci¬ 
ton wave function confined in a single layer. Our results 
exemplify the control of spatial dimensionality of the val¬ 
ley excitations in 2H- and 3R-MA'2, which may help to 
develop new valleytronics devices. 


IV. DISCUSSIONS 

The present results establish the 2D nature of the val¬ 
ley exciton in bulk M 0 S 2 with the 3R stacking. Note 
that an apparently similar confinement effect has also 
been discussed for the 2H stacking as the spin-layer lock¬ 
ing effect mm- Since the magnitude of confinement 
is governed by the factor X/^/X 2 +t% with A being the 
spin-orbit coupling coefficient [3T], its strength essentially 
relies on atomic number of M. Distinctively, the layer- 
confinement effect in the 3R stacking is a quantum inter¬ 
ference effect caused by the symmetry of the crystal and 
therefore is guaranteed to be relevant in general MX 2 
systems with the same structure. The complete confine¬ 
ment of the 3R valley states should also provide us with 
an arena of exploring novel interplay of low dimensional¬ 
ity and valley-physics. 

Compared to the 2H or monolayer MoS 2 , the 3R poly¬ 
type has a significant merit that valley polarized carriers 
are more robust in the latter form. The carriers are pro¬ 
tected from relaxation processes resulting from interlayer 
transitions, contrary to the 2H systems. In the previous 
work 133] . we have found an indication that interlayer 
transitions are much suppressed and the valley excitons 
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Appendix A: FITTING OF THE IMAGINARY 
PART OF THE DIELECTRIC FUNCTION 

The imaginary part of the dielectric function deduced 
from the measured reflectivity is fitted by a sum of 
Lorentzian peaks. The fits shown in Fig. [5] describe the 
experimental spectra well. The contributions from the 
higher energy transitions are also included. 
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FIG. 5: (a)-(b) Fit of the imaginary part of the dielectric function in 2H- and 3 R-M 0 S 2 , respectively. The dielectric function 
deduced from the reflectivity (black lines) can be fitted by a sum of Lorentzian peaks (red lines). The green lines show the 
contribution of the respective peaks. 
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